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We calculate nonequilibrium quasiparticle and phonon distributions for a number of widely- 
used low transition temperature tlrin-film superconductors under constant, uniform illumination 
by sub-gap probe and pair-breaking signal photons simultaneously. From these distributions we 
calculate material-characteristic parameters that allow rapid evaluation of an effective quasiparticle 
temperature using a simple analytical expression, for all materials studied (Mo, Al, Ta, Nb, and NbN) 
for all photon energies. We also explore the temperature and energy-dependence of the low-energy 
quasiparticle generation efficiency g by pair-breaking signal photons finding g fts 0.6 in the limit of 
thick films at low bath temperatures that is material-independent. Taking the energy distribution of 
excess quasiparticles into account, we find g —> 1 as the bath temperature approaches the transition 
temperature in agreement with the assumption of the two-temperature model of the nonequilibrium 
response that is appropriate in that regime. The behaviour of g with signal frequency scaled by the 
superconducting energy gap is also shown to be material-independent, and is in qualitative agreement 
with recent experimental results. An enhancement of g in the presence of sub-gap (probe) photons is 
shown to be most significant at signal frequencies near the superconducting gap frequency and arises 
due to multiple photon absorption events that increase the average energy of excess quasiparticles 
above that in the absence of the probe. 

PACS numbers: 74.40.Gh, 74.78.-w, 29.40.-n, 74.25.N- 


I. INTRODUCTION 

The energy distribution of the quasiparticle excita¬ 
tions of a superconductor determines its electricaP and 
thermal transport properties!- 2 A theoretical understand¬ 
ing of devices including kinetic inductance detectors 
(KIDs)pHS superconducting tunnel junctions^ supercon¬ 
ducting qubits?® SQUID based parametric amplifiers 
and mixersp^ and quantum capacitance detectors 11 re¬ 
quires a description of the quasiparticle system response. 
Therefore, a method for calculating the excited, possibly 
nonequilibrium distribution of quasiparticles is central to 
any device model. 

Nonequilibrium superconducting detectors rely on pair¬ 
breaking to create their detected signal. We are partic¬ 
ularly interested in thin-film superconductors fabricated 
on a dielectric substrate that is cooled to a sufficiently 
low temperature (the bath temperature) Tj,. Typically 
Xb ~ 0.1 T c where T c is the superconducting transition 
temperature. An interacting signal photon of energy 
/^signal, where h is Planck’s constant and Signal the sig¬ 
nal frequency, directly breaks a condensate pair in the 
superconductor provided hv s i gna i > 2A where A is the 
superconducting energy gap. Due to the high density of 
states close to the gap, and for moderate energy photons, 
the interaction creates a distribution of excited quasiparti¬ 
cles with peakJh 2: at E = A and E = hv. s i gna i — A where E 
is the quasiparticle energy. E is related t o the underlying 
Bloch state energy e by E = \J e 1 + A 2 . Higher energy 
signal photons will directly release an atomic electron via 
the photoelectric effect. In this case effects associated 
with localized heating and gap-reduction or “hotspots” 


may become important. In this work we concentrate 
on modelling the detection of moderate energy photons 
hz4i g nai < where these effects are insignificant. The 
photon energies considered are particularly important 
to understand the responsivity, sensitivity and signal- 
to-noise of currently-deployed and planned astronomical 
instruments performing measurements in the frequency 
window of 0.1 to 10 THz. In the case of KIDs, the nonequi¬ 
librium state created by the interaction is monitored using 
a probe of energy hv plo b e <C 2A, where t'probe is the probe 
frequency, so an additional drive term needs to be included 
in the detailed analysis. 

An important consideration in any nonequilibrium su¬ 
perconducting detector operating at low reduced tem¬ 
peratures Tb/T c <C 1 is that the detected signal is most 
influenced by the presence of low-energy quasiparticles 
because the relaxation of the primary excitation occurs 
on time-scales that are much shorter than the ultimate 
relaxation time of the low-energy excess. (Here we define 
low-energy quasiparticles to have E < 3A.) We assume 
the signal photons interact with 100% efficiency i.e. we 
ignore for example any optical coupling efficiencies. A pri¬ 
mary excitation with E = liv s i gn ai — A relaxes towards the 
gap emitting a phonon. If this emitted phonon has energy 
fl > 2A, an additional pair may be broken enhancing the 
number of low-energy excess quasiparticles and thus the 
detected signal. But in a thin-film this emitted phonon 
may be lost into the substrate, reducing the number of 
low-energy excitations. The average number m of low- 
energy quasiparticles resulting from the interaction of a 
single high-energy photon can be quantified in terms of a 
quasiparticle generation efficiency (or quasiparticle yield) 
g. If the low-energy quasiparticles have average energy 
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(E qp ) then the efficiency 77 = m(E qp )/his s i gna i. Often it is 
assumed (E qp ) = A . 13 

Phonon loss means 77 < 1. Most existing work assumes 
a value of 77 ss 0.6 for all materials, photon frequencies and 
bath temperatures]^- At high temperatures T ss T c our 
description of the energy-cascade no longer applied ® 15 
because thermal phonon energies become comparable with 
A, and scattering and recombination occur on comparable 
timescales. Then a two-temperature model describing the 
quasiparticle and phonon systems that assumes 77 = 1 is 
most appropriate. Here we show that the full nonequilib¬ 
rium calculation yields 77 —> 1 as Tf,/T c —► 1 in agreement 
with the assumption of the two-temperature model. 

Previous work modelled the effect of incoming energy 
at low temperatures with an effective temperature 16 or 
as an effective chemical potential 17 ! for the quasiparti¬ 
cle distribution. To go beyond this analysis, Chang and 
Scalapino 18 19 derived a set of coupled nonlinear kinetic 
equations describing the interactions of quasiparticles and 
phonons to find their respective energy distributions f{E) 
and n(fl). They solved for the nonequilibrium distribu¬ 
tions resulting from various drive terms, including photon 
and phonon injection. A number of investigations have 
explored the effect of very high energy pho tons (x-ray or 
optical photons) on infinite superconductors 20 ® and thin 
filmJ 22 ! calculating a quasiparticle generation efficiency 
of 77 ~ 0.6, ignoring Q > 2A phonon loss. Another ap- 
proaclfBO considered the time-evolution following local 
energy deposition into the thin film. Kozorezov et al J 23; 
considered the energy downconversion process after ab¬ 
sorption of a high energy photon for a variety of materials 
and concluded that the materials can be categorized into 
three different classes, with the low energy-gap supercon¬ 
ductors all having 77 ss 0 . 6 . 

We have previously reported® a numerical approach 
that solves the coupled kinetic equations describing the 
quasiparticle and phonon distributions f(E) and n(fl) for 
a superconducting thin-film driven by a sub-gap probe!® 
and including the effect of an additional pair-breaking 
signal , 25 for Al thin-films at low temperatures T ~ 0.1 T c . 
We have also reported detailed calculations of KID charac¬ 
teristics (quality factors and quasiparticle lifetimes as func¬ 
tions of readout power and Xj,) that were compared to pre¬ 
cise measurements of Al resonator behaviour finding good 
agreement.® More recent measurements on KIDs 27 of the 
quasiparticle number and recombination time-dependence 
on readout power appear to show that they must be ex¬ 
plained using nonequilibrium quasiparticle distributions. 
Other recent measurements with a Ta KIE® report a de¬ 
tector response that has the same qualitative features as 
the energy dependence of 77 for Al described in Guruswamy 
et al. 25 Here we apply the method to a number of other 
technologically important low temperature superconduc¬ 
tors and discuss the physics involved. We concentrate 
in particular on the refractory elements and extend the 
earlier analysis to bath temperatures approaching the 
superconducting critical temperature T c . We study Mo, 
Al, Ta (in its higher T c ~ 4.4K form), Nb and NbN. In 


each instance we use measured material properties that 
characterize films that we deposit by sputtering under 
ultra-high-vacuum. Our results are applicable to many 
different device designs, as only the phonon escape time 
to the substrate 77 is device-geometry dependent. The 
method assumes that phonon pair-breaking, quasiparticle 
recombination, and electron-phonon scattering are the 
only significant interaction processes. This assumption is 
investigated in section II B| for typical thin-films intended 
for detector applications. Section[H]reviews superconduct¬ 
ing parameters for the materials discussed, considers the 
relative contribution from electron-electron scattering for 
the different materials, describes the numerical method, 
and also describes how 77 is calculated for a pair-breaking 
signal. Section m presents an analytical superconductor 
cooling model that can be used to calculate the effec¬ 
tive quasiparticle temperature T*n for both sub-gap and 
pair-breaking photons. Section [TV] describes detailed nu¬ 
merical results for both sub-gap and pair-breaking photon 
interactions including the effect of the effect of changes 
in the bath temperature. Section |V| summarizes the work 
done. 


II. METHODS 


A. Superconducting parameters 


The superconducting parameters needed for modelling 
are the characteristic quasiparticle lifetime 


Zi(0)7i 

r °" 27r6(fc B r c )3 ’ 

and the characteristic phonon lifetime 


hNj on 

47r 2 iV(0) bs (a 2 )A o ’ 


(1) 

( 2 ) 


as defined by Kaplan et al.® IV(0) and lV(0)b s are the 
mass-enhanced and bare single-spin density of states 
at the Fermi energy ep respectively, where N(0) = 
Zi(0)N(0)b s , Won is the ion density, and A 0 the zero 
temperature superconducting energy gap. Z\( 0) is the 
electron-phonon renormalization factor Z\( 0) = 1 + A, 
where A = 2 f Q dfla 2 (fl)F(fl)/fl is the dimensionless 
electron-phonon coupling strength. a 2 (fl)F(fl) is the 
Eliashberg function, a 2 (fl) is the electron-phonon inter¬ 
action matrix element and F(Q) is the phonon density of 
states. The material-dependent parameter b is calculated 
within a Debye model such that bit 2 = a 2 (il)F(fl) for low 
phonon energies D ~ 2 Aq. The averaged value of the in¬ 
teraction matrix element is (a 2 ) = 1/3 f Q dfl a 2 (it)F(fl). 
Values of b and (a 2 ) for Al, Ta and Nb are already tabu¬ 
lated.^ For Mo we use the point-contact measurements 
of Caro et al.® while for NbN we use the tunnelling data 
of Kihlstrom et al. 31j and the transport property mea¬ 
surements of Semenov et al.® We note that the precision 
of our calculations of b are dependent on the accuracy 
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within which we can interpret the relevant low-frequency 
parts of the reported a 2 (£l)F(Q) data. 


Material 

Mo 

Al 

Ta 

Nb 

NbN 

A 0 (peV) 

140 

180 

700 

1470 

2 560 

u s&p (GHz) 

68 

87 

339 

711 

1240 

T c /K 

0.92 

1.18 

4.6 

9.67 

16.8 

b (10 -4 meV -2 ) 

2.28 

3.17 

17.3 

40 

47 

(a 2 ) (meV) 

1.62 

1.93 

1.38 

4.6 

4.99 

A 

0.42 

0.43 

0.69 

1.84 

1.46 

to (ns) 

1310 

438 

1.78 

0.15 

0.02 

To (PS) 

231 

260 

22.7 

4.17 

5.98 

Tec/ Tq 

1.04 

1.04 

0.91 

0.94 

1.01 


Table I. Characteristic quasiparticle and phonon lifetimes, 
and associated parameters. Data from Gladstone et aZ.^“, 
Kaplan et al) or Zehnder 13 unless otherwise specified. For 
Mo and NbN we use measurements of a 2 (Q)F( fi) from Caro 
et aZ.®® and Kihlstrom et al. 31 ^ respectively. t£ c is the value 
of t$ required for energy conservation calculated from (JsJ) . 
v g a p = 2 Aq /h is the gap frequency. 


The characteristic times calculated are listed in ta¬ 
ble H| We assume the weak-coupling relationship Ao = 
1.76 k B T c for all materials to determine T c in table [TJ The 
assumption is reasonable for Mo Al and Ta, but gives a 
slightly higher T c compared to experiment for Nb (typ¬ 
ically T c ~ 9.3 K) and for NbN (T c ~ 15 K). However 
we do not expect the assumption to significantly affect 
the physics of the conclusions drawn; the differences in¬ 
troduced by strong coupling are reduced becau se of the 
short lifetimes of states with energy E A. 18 2 ^r^. is the 
value of Tq required for energy conservation calculated 


from and discussed in more detail in section IIC 


B. Electron-electron vs electron-phonon scattering 


Insight into the contribution from electron-electron e—e 
scattering in the down-conversion process is most-easily 
found by considering the relavant scattering rates T c . For 
the electron energies of interest e <C fZp, we assume that 
electron-phonon (e—p) inelastic scattering dominates over 
(e — e) scattering. We previously investigated 2 ® to what 
extent this assumption is valid for low-resistivity Al films. 
Now we consider the other low-T c superconductors using 
material parameters typical of those that we measure 
experimentally for the thin-films we deposit by magnetron- 
sputtering in ultra-high vacuum. 

In the normal-state the e — p scattering rate t~_/ p i^U 

r e 

r-_ 1 p (e) = -j clfl« 2 (n)F(9). (3) 

Assuming a 2 (fZ) is constant at low energies, a Debye 


model for phonon^Hl leads to 


Te—p(^) 


e 


3 


3 T 0 (k B T c ) 3 


3T 0 (k B T c ) 3 


if e < Dp 
if e > Dp 


( 4 ) 


The electron-electron scattering rate in clean films is 
estimated with the Landau-Pomeranchuk formula 2 ® 


-i 


-e(e) = 


2 1/2 
t 'S 

he^Y. 96’ 


( 5 ) 


where r s is the radius containing one electron charge 
divided by the Bohr radius (approximated as 1). For thin 
resistive films the e — e rate can be significantly enhanced 
so that 8 ® 


T e— e( e i Rsq) 


e 2 R sq e _! for 


( 6 ) 


where R sq = p/d is the sheet resistance of the thin-film, 
p is the resistivity and d is the thickness.^® These normal- 
state calculations significantly overestimate the e — e scat¬ 
tering rates in the superconducting state at low tempera¬ 
tures Tf, T C W^ Figure [l] plots the normal-state scatter¬ 
ing rates for a 40 nm thick Mo Him and we have used the 
thin-film resistive e — e scattering rate ©■ We identify 



Figure 1. Energy dependence of normal state scattering rates 
in Mo. Electron-phonon rate r~_ p (e) (solid green), clean 
limit electron-electron rate T~_ e (e) (dotted blue), and thin- 
film resistive T~_ e {e, R aq ) (dashed red). Calculations are for a 
40 nm Mo film with p = 9.2 x 10~ 8 fim. 

the following characteristic energies for each material: cf', 
Qd; and a low-energy crossover E c , below which the e — e 
scattering rate exceeds the e—p rate. E c j A is highest 
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Material 

I— 1 

0 

1 

00 

£) 

d (nm) 

e F (eV) n D (: 

meV) E c / A 0 

Mo 

9.2 

40 

9.32 

39.6 

6.1 

aP 

0.8 

35 

11.6 

36.2 

1.42 

Ta (bcc, epi) 

4.1 

100 

9.5 

20.7 

0.23 

Nb 

OO 

bo 

100 

6.18 

23.7 

0.14 

Nt Dj39140| 

250 

100 

15.6 

28.4 

0.3 


Table II. Characteristic material energies related to scattering. 
E c is the energy at which the electron-phonon scattering 
rate becomes greater than the thin-film resistive electron- 
electron scattering rate. Data from Gladstone et al .p Kaplan 
et al . 29 or Kozorezov et aZ.^® unless otherwise referenced. 
The material resistivities and thicknesses used are measured 
values from thin-films deposited by our group. All data is for 
polycrystaline films except for Ta where we assume parameters 
typical of the higher- T c bcc form. 


for Mo for the material parameters chosen, as shown in 
table HIl 

During the energy down-conversion, the inelastic scat¬ 
tering rate is only relevant above energies E > 3A where 
additional pair breaking is possible. We find that E c is 
below 3A for all materials except Mo. Figure [l] shows 
that e — p scattering dominates in that case for e ss 6 A to 
10 4 A for the modelled Mo film, and e — e scattering may 
contribute above and below this energy range. However, 
the e — e scattering rate is considerably suppressed in the 
superconducting state at low T/T C W -1 We conclude that 
there is negligible contribution from e — e scattering for 
most of these materials, thicknesses, temperatures and 
signal photon energies. For Mo we note that additional 
e — e scattering may increase the quasiparticle generation 
efficiency over the results presented here by providing an 
additional pair-breaking mechanism especially near T c . 


C. Numerical method 


We solved the coupled kinetic equation^ 1 ®! describing 
the quasiparticle and phonon distributions numerically to 
find the steady-state driven distributions, f(E) n( f2) for 
the quasiparticles and phonons respectively, using Newton- 
Raphson iteration. Drive terms for pair-breaking and sub¬ 
gap photoni 4l -lwere included as necessary and the numeric 
pre-factors of these terms scaled to match the absorbed 
powers, Pprobe for sub-gap photons and P s i gn ai for the 
pair-breaking signal. The numerical method was iterated 
until the solutions conserved energy to within £ = 10 -5 . 

This was quantified as £ = + £|_ b , where £ 9P _0 is 

the relative difference between the photon power absorbed 
and power flow from quasiparticles to phonons, and £;,p-b is 
the relative difference between the photon power absorbed 
and the power flow from phonons to substrate (held at 
a temperature T b ). The nu merical met hod is outlined in 
more detail in Goldie and WithingtonP 


For each material, the temperature dependent energy 


gap A(T) was calculated using the BCS gap equation, 

i _ i-2/(p,t) m 

A^( 0 )Vbcs Ja(t) ^E 2 - A(T) 2 ’ 

where f(E, T) = (exp (P/fc^T) + 1 ) _1 is the Fermi-Dirac 
distribution, and 7V(0)Vbcs is the dimensionless electron- 
phonon coupling constant, calculated from the zero tem¬ 
perature form of 0 as -ZV(0)Vbcs = 1/sinh - 1 (fl£>/Ao). 

The numerical method converges on a solution that 
minimizes £ provided the material parameters satisfy 


27riV(0)T 0 ^Aof2®, 

9IWo(W 


a consequence of matching the low energy part of the mea¬ 
sured Eliashberg function to a Debye-model approxima¬ 
tion that determines b. In our previous work we ensured 
that T c satisfied 0 pll but then the T c obtained did not 
precisely agree with 0 (although the difference is small 
< 1%). In this work, since we are also interested in the 
temperature dependence of the response we calculated T c 
from Q and defined an energy-conserving characteristic 
phonon lifetime t£ c to satisfy (18]) . Given the difficulty of 
interpreting the low-energy a 2 (T l)F(Q) data that we and 
indeed Kaplan et al .P and others have noted we consider 
this approach entirely reasonable. The ratios t£ c /tq are 
given in table [TJ Satisfyingly the ratios are within 10 % 
of unity for all materials studied. In our calculation of b 
for Mo, our error would be of this order given the diffi¬ 
culty of determining the low-energy a 2 {p,)F(fl) from the 
available data. We note that other estimates of the char¬ 
acteristic times exist (see for example Parlato et al .^) 
although those estimates do not necessarily satisfy the 
energy-conserving requirement. 

In the numerical code the energy-bin distribution for 
both the quasiparticles and phonons was also changed 
when calculating bath temperature dependencies. As T b 
is increased, the phonon-bath power flow P(f 2 ) 0 _ b above 
fl = 2A increasingly resembles the difference between 
two thermal distributions, and an increasing fraction of 
power is carried by high energy phonons, as shown in 
figure [2] Therefore at high temperatures, the phonon 
distribution must be accurately calculated to these higher 
energies. The energy of the last phonon bin fl max signif¬ 
icantly affects the energy conservation of the solutions. 
Figure [ 2 ] (inset) shows that to have a relative power flow 
difference £^-6 < 10~ 5 , required fl max > 15 fesTj, for 
T b /T c = 0.95. 

However, the bin width could not be made too large 
as the energy gap A (T), readout photon energy hv p , and 
signal photon energy hv s had to be rounded to the nearest 
multiple of the bin width, such that the photon-induced 
peaks in the distributions occurred within well-defined en¬ 
ergy bins. We chose to model the phonons with N = 2500 
uniformly-sized bins with maximum phonon energy S 2 max , 
and similarly the quasiparticles with maximum energy 
Fmax = Umax + A (T). Bin width was therefore variable, 

















5 



tiond^ for T $ and r ?p results iiP^l 


V ~ (E qp ) 


WLial -^ 2 


probe/ 


signal 


2 R 

r +ja ’ 


(9) 


where R and ft are the distribution-averaged recombina¬ 
tion rate and pair-breaking rates respectively. We define 
(E qp ) the average energy of the excess low-energy quasi¬ 
particles where 


, _ f~Ep(E)(f(E)-f(E,T))dE 

QP) f~p(E)(f(E)-f(E,T))dE ' 


( 10 ) 


At low temperatures (E qp ) « A as usually assumed: tak¬ 
ing account of the detailed energy distributions allows i) 
to be determined at arbitrary temperatures. 


III. SUPERCONDUCTOR COOLING MODEL 


Figure 2. Plionon-batli power flow P (at Tb = 0.1 T c 
(solid red) and Tb = 0.95 T c (dashed blue). Both for NbN, 
Pprobe — 2 X 10 W m , Psignal — 0, TiJTq — 1, and hv probe ~ 
16 fieV. Inset of relative phonon-bath power flow difference 
against the energy of the maximum bin flmax/fcsT),, at 
the same bath temperatures Tb = 0.1 T c (red x) and Ti, = 
0.95 T c (blue +). 


We used the steady-state driven distributions that were 
calculated to determine parameters for the thin-hlm su- 
percond uctor cooling model outlined in Goldie and With-| 


ington [23 This model is an analytic expression relating 


Tj( r to the power flow P between the quasiparticle and 
phonon systems (the latter assumed to be the substrate 
phonons) for a given Tb <C T c and material. 


depending on material and temperature. The quasiparti¬ 
cle density of states broadening parameter 7 was recalcu¬ 
lated for each simulation, changing for different material 
and bin width combinations. It was chosen to minimize 
the difference between the quasiparticle number calcu¬ 
lated as the integral 41V(0) dE Pbcs(E, A)f(E, T b ), 
where pbcs(E,A) = E/y/E 2 — A 2 , and the equivalent 
sum over the discretized distribution with the broadened 
density of states calculated with E —► E + i'y. In all cases 
the approach worked well. 


D. Low-energy quasiparticle generation efficiency 
for a pair-breaking signal 

For a single pair-breaking photon we have defined a low- 
energy generation efficiency r] = m(E qp )/hu s i gna i. Here 
we describe the calculation of p for a constant input 
signal power P s i gn ai in a thin-film that may also be driven 
by a non-direct pair-breaking constant readout power 
Pprobe (such as used for KID readout). The absorbed 
power creates a number of excess quasiparticles fV excess = 
A^ignai — lVp ro be such that the quasiparticle system total 
energy has increased by E excess = E si gnal - P pro be, where 
A^ignai and P s i g nai are the quantities for distributions 
driven by the signal and probe together. In the steady 
state, photons absorbed at a rate T$ create low-energy 
quasiparticles at a rate r ?p = mr$. 

Solving a modified set of Rothwarf-Taylor rate equa- 


P = 


p{P : u) 1 + Tl/Tpb ' 


^ f -2A {T* n ) \ rr 

r " exp (^r) 1 p 


-2A(r t ) 

ksTb 


• ( 11 ) 


E s is a material-dependent constant, r p b is the phonon 
pair breaking time, ~ Tq for Tb <C T c . T* n is the effec¬ 
tive temperature calculated to characterize the driven 
nonequilibrium quasiparticle distributions so that N qp = 
47V(0) /^° d E p(E)f(E, Tjy). 77 (P, v) depends on the drive 
(probe or pair-breaking signal). Goldie and Withington^ 
showed for a sub-gap probe rj(P, t'probe) = V 2 A the frac¬ 
tion of phonon-bath power flow carried by excess phonons 
with energy f 1 > 2A(T b ) that depends on P, t'probe and 
Ti. Here we show that (11) can also be applied for di¬ 
rect pair-breaking when 77 (P, Signal) = 77 as defined in 
section III D1 


IV. RESULTS 


A. Effects of a sub-gap probe 


Here we show and discuss results for modelling of the 
effect of a sub -gap probe. The data-points (crosses) 
in figure |3(a) show the calculated T^ from the full 
non-equilibrium solution as a function of P pro be- The 
solid line shows T ^ evaluated with © and from this 





































6 





Figure 3. |(a)| Effective quasiparticle temperatures T% as a 
function of P pro b e ; full nonequilibrium calculation (blue +), 
and a fit using (111 (red solid line), (b) Fraction of phonon- 
bath power flow carried by fi > 2 A excess phonons as a 
function of P pro b e ; full nonequ ilibr ium calculation (blue +), 
and piecewise fit (red line) to (12 I. The dotted vertical line 
marks Po. 


Figure 4. Parameters for piecewise fits to 772 A for Nb (blue 
x), Ta (green □), and A1 (red O) calculated at a range of 
frequencies hv pro b e - (a) Pq against scaled frequency, (b) 770 


(constant section) against scaled frequency, (c) 771 (slope for 
linear section) against scaled frequency. 


Ml Distribution-averaged phonon pair-breaking emphasize the commonality. 


time T p b- The characteristic phonon lifetime Tq is marked 


with the dotted line, (d) 772 A as a function of phonon escape 


Material A 0 (/reV) E a (Wm^KT 1 ) P 0 (Wm“ 3 ) 


2 x 10 3 Wm 3 , hv probe = 16/reV, and Ti/rfi = 

l. 


Mo 

140 

1.42 • 

10 10 

2.96- 

10“ 3 




A1 

180 

3.23- 

10 10 

2.83- 

10“ 3 




Ta 

700 

1.17- 

10 14 

1.30- 

10 1 

we also calculate £,. The calculations 

are for 

Nb, 

Nb 

1470 

3.83- 

10 15 

3.02- 

10 3 

with P probe — 2 X 10 W ID. , /pprobe — 

16 /reV, 

and 

NbN 

2 560 

1.29 • 

10 16 

3.17- 

10 4 


ti/ t q = 1. The analytical model is an excellent approx¬ 
imation to the full nonequilibrium calculation provided 
the power-dependence of 772 a shown in figur e |3(b)| is taken 
into account. The solid line in figure [3~(b)| is a piecewise 
fit such that 


Table III. Superconductor cooling model material parameter 
E s , and knee power Po at hv pro be = 16 /reV and ti /tq = 
1, obtained from fits to results of the full nonequilibrium 
calculation. 


/px jvo ifP<P 0 , 

??2A (P)=< wp/ px, -rp^p ( 12 ) 

771 ln(P/P 0 ) + Vo if P > Po- 


Po characterizes a “knee” probe power below which 772 A 
is constant and 771 characterizes the energy dependence of 
772 A at highe r probe powers. Both are material dependent. 
Figure 3(c) shows that r p i, ~ Tq is a good approximation 
over the range of interest. Figure [3(d)| shows that 772 a 
has no significant dependence on the phonon escape time 
p. 

The dependence of the fit parameters Po, r/ 0 , and rji on 
probe photon energy is shown in figure [4] This includes 
calculations for Nb (blue x), Ta (green □), and A1 (red 
O)- Scaling the energy of the photon by A 0 demonstrates 
the monotonic behaviour for all the materials. Pq shown 


in figure 4(a) is scaled by the material-dependent T, s to 


Table [III] summarizes parameters derived from the mod¬ 
elling to describe sub-gap photon interactions in all of the 
materials studied. The table shows £ s and knee param¬ 
eter Po at hv probe = 16 77 eV and ti/tq = 1. The values 
of £ s shown can also be used for estimates of for a 
direct pair-breaking signal as we will discuss in the next 
section with an appropriate energy-dependent 77. 


B. Effects of direct pair-breaking 

Here we discuss results for a direct pair-breaking sig¬ 
nal with 2A < hv s [g na \ < 10A. Figure [5] shows how 77 
varies with signal photon energy. 77 decreases from unity 
in the range hv s i gna i = 2A to 4A, as phonons emitted 
by quasiparticles scattering towards the gap escape the 
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Figure 5. Quasiparticle generation efficiency r/ po wer against sig¬ 
nal photon energy hv a i gns i for different materials. Calculated 
for T b = 0.1 T c , Pp ro be = 2 x 10 J Wm -3 , hv pro be = 16/reV, 
Psignai = 2Wm' 3 , and ti/t$ = 1. 


thin film. Above hv s i gn ai = 4A, the phonons emitted in 
scattering can break additional Cooper pairs, increasing 
the quasiparticle generation efficiency. The magnitude 
of the increase depends on the phonon trapping factor ! 25 
Measurements of the response of a Ta KID--® at signal 
frequencies close to the gap energy show these character¬ 
istic features. The figure shows that when hv signal, Tb 
and T[ are scaled by the relevant material parameters, 
the quasiparticle generation efficiency is the same for all 
materials. Kozorezov et al) 23 predicted that the quasi¬ 
particle generation efficiency could be categorized into 
three classes of superconductors. The superconductors 
included in this work, which fall within either the second 
(Ta, Nb, NbN) or third (Al, Mo) class, have a material- 
independent r) ss 0.6. This conclusion differs somewhat 
from that of IZehnderP ^ We understand this as due to 
the current work examining the steady-state response to 
constant incoming power, and scaling all relevant param¬ 
eters by material characteristics, whereas |Zehnderf s work 
examined the time-dependent response including quasi¬ 
particle out-diffusion from a localized hot-spot, and chose 
a fixed cutoff time of 10 ns to calculate the quasiparticle 
generation efficiency unsealed by material. We emphasize 
that for the low energy photons considered in the present 
work localized hot-spot formation (and consequent gap 
suppression) does not occur for the signal powers studied. 

The other effect shown in figure [5] is the enhancement of 
77 by the readout power due to an increase in the average 
energy of generated quasiparticles. For hv s i gna i — 2A, it 
can be seen that 77 > 1 which may seem unreasonable. 
The enhancement is due to multiple photon absorption 
by the quasiparticles i.e. both signal and probe, and 


the effect is most-pronounced near hv s i gn ai = 2A. This 
results in a higher average excess quasiparticle energy 
when the readout power is present alongside the signal 
compared to when only the signal is present, and also 
produces the small variation in 77 between materials. The 
absorbed readout power P pro be and readout frequency 
hv probe ar e not scaled by material parameters in this 
figure, unlike the signal. When calculated with P pro be = 0, 
the signal quasiparticle generation efficiency is identical 
for all materials and is unity at hv sig na i = 2A. 



T b /T c 


Figure 6 . The quasiparticle generation efficiency 77 as a func¬ 
tion of reduced temperature T b /T c . The calculation is for for 

Al, with Pprobe = 0 , Psignal = 2 Wm“ 3 , hi 'signal = 10A(T&), 
and two values of the phonon escape time. The inset shows Ty 
calculated using the full nonequilibrium model (red O) an d 
using <[TT|) (blue line) as a function of P s i g nai with T b = 0.1 T c 
and tz/tq = 1 . 


Figure [ 6 ] shows the variation in 77 as the bath tempera¬ 
ture is changed. As T b increases, the scattering between 
the more numerous thermal quasiparticles and phonons 
determines the structure of the driven distribution, rather 
than scattering between the excess quasiparticles and 
phonons. This means the power flow from phonons to 
bath does not have the same peaked structure and instead 
more closely resembles the difference between two thermal 
distributions, as earlier shown in figure [2j This type of 
excess quasiparticle distribution means the average energy 
of generated quasiparticles increases without increasing 
the recombination and pair breaking rates. The figure 
also indicates temperatures for which ksTb = 2A(T b ) and 
= A (T b ), above which we expect thermal and nonequilib¬ 
rium distributions to interact strongly. As T b /T c —► 1 we 
see that 77 —1 1 in agreement with the two-temperature 
model valid in that regime. We find that T^ is well- 
accounted for by ( 11 ), with the value of 77 shown, for 
Tb/T c <0.7. The inset shows T ^ calculated using 


















V. CONCLUSIONS 


for ti/tq = 1 as a function of -P s i gn ai with Tb = 0.1 T c . 



Figure 7. Quasiparticle generation efficiency r/ for materials 
studied as a function of phonon escape time ti/tq . Calculation 
uses Tb = 0.1 T c i -fprobe — 0, -/“signal — 2 W m 3 , and /iz/ s ig na l = 

10 A. 


Figure [7] shows the dependence of rj on the phonon 
trapping factor r;/r p b for Tb = 0.1 T c . As the phonon 
escape time ti/tq increases, the quasiparticle generation 
efficiency rj —> 0.6 for all materials examined. This con¬ 
clusion is in agreement with earlier Monte Carlo calcula¬ 
tions 20 22 for infinite superconductors or where Q > 2A 
phonon loss was ignored. 


We have described solutions of the coupled kinetic equa¬ 
tions that calculate steady-state nonequilibrium quasipar¬ 
ticle and phonon distributions in a number of techno¬ 
logically important superconducting thin-films driven by 
sub-gap and pair-breaking photons. In particular we con¬ 
sider low energy signal photon interactions hi/ s - lgna \ < 10A, 
where localized heating and gap-suppression can be ig¬ 
nored. We have calculated numerical parameters for these 
superconductors that can be used in a simple analytical 
expression to describe the power-flow between quasiparti¬ 
cles and phonons. This expression allows straight-forward 
estimates of the effective quasiparticle temperatures for 
both sub-gap and pair-breaking drives to approximate 
the nonequilibrium behaviour without resorting to a full 
numeric solution of the coupled kinetic equations. The an¬ 
alytical expression is shown to give a good account of full 
solutions of the detailed equations for a wide range of pow¬ 
ers and bath temperatures. This is relevant for predicting 
the behaviour of thin-film nonequilibrium superconduct¬ 
ing detectors. We defined a low-energy quasiparticle 
generation efficiency for constant absorbed pair-breaking 
power and calculated detailed numerical solutions as a 
function of photon energies and bath temperature. We 
have shown that key parameters determining quasipar¬ 
ticle generation efficiency are the signal frequency and 
the phonon trapping factor. The enhancement of signal 
quasiparticle generation efficiency by absorbed readout 
power is demonstrated and explained as multiple-photon 
absorption (signal and probe) by the quasiparticles and 
Cooper pairs. We estimate this effect can result in an 
increase r] of up to 20% at signal photon energies of 
^signal = 2A. The most sensitive detectors currently 
in development may also be able to distinguish between 
phonon trapping factors using results presented here. For 
the low-energy photon interactions studied we also show 
that the quasiparticle generation efficiency is material 
independent. 
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